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Abstract 

In quantum gravity, several approaches have been proposed until now for the quan¬ 
tum description of discrete geometries. These theoretical frameworks include loop 
quantum gravity, causal dynamical triangulations, causal sets, quantum graphity, and 
energetic spin networks. Most of these approaches describe discrete spaces as homo¬ 
geneous network manifolds. Here we define Complex Quantum Network Manifolds 
(CQNM) describing the evolution of quantum network states, and constructed from 
growing simplicial complexes of dimension d. We show that in d = 2 CQNM are ho¬ 
mogeneous networks while for d > 2 they are scale-free i.e. they are characterized by 
large inhomogeneities of degrees like most complex networks. Prom the self-organized 
evolution of CQNM quantum statistics emerge spontaneously. Here we define the 
generalized degrees associated with the d-faces of the d-dimensional CQNMs, and we 
show that the statistics of these generalized degrees can either follow Fermi-Dirac, 
Boltzmann or Bose-Einstein distributions depending on the dimension of the d-faces. 
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Several theoretical approaches have been proposed in quantum gravity for the description 
and characterization of quantum discrete spaces including loop quantum gravity M, causal 
dynamical triangulations HIS], causal sets PE], quantum graphity IHIDI, energetic spin 
networks mm, and diffusion processes on such quantum geometries PI- In most of these 
approaches, the discrete spaces are network manifolds with homogeneous degree distribution 
and do not have common features with complex networks describing complex systems such 
as the brain or the biological networks in the cell. Nevertheless it has been discussed [H] that 
a consistent theory of quantum cosmology could also be a theory of self-organization [TSl EB] , 
sharing some of its dynamical properties with complex systems and biological evolution. 

In the last decades, the held of network theory [IMI] has made signihcant advances in 
the understanding of the underlying network topology of complex systems as diverse as the 
biological networks in the cell, the brain networks, or the Internet. Therefore an increasing 
interest is addressed to the study of quantum gravity from the information theory and 
complex network perspective [221123] • 

In network theory it has been found that scale-free networks [23] characterizing highly 
inhomogeneous network structures are ubiquitous and characterize biological, technological 
and social systems [IMQ]. Scale-free networks have hnite average degree but inhnite huc- 
tuation of the degree distribution and in these structures nodes (also called ’’hubs”) with 
a number of connections much bigger than the average degree emerge. Scale-free networks 
are known to be robust to random perturbation and there is a signihcant interplay between 
structure and dynamics, since critical phenomena such as in the Ising model, synchronization 
or epidemic spreading change their phase diagram when dehned on them [25] l26] . 

Interestingly, it has been shown that such networks, when they are evolving by a dynamics 
inspired by biological evolution, can be described by the Bose-Einstein statistics, and they 
might undergo a Bose-Einstein condensation in which a node is linked to a hnite fraction of 
all the nodes of the network [27]. Similarly evolving Cayley trees have been shown to follow 
a Fermi-Dirac distribution [28] l2^. 

Recently, in the held of complex networks increasing attention is devoted to the charac¬ 
terization of the geometry of complex networks [SOHSO]. In this context, special attention 
has been addressed to simplicial complexes [3QH33], i-e. structures formed by gluing together 
simplices such as triangles, tetrahedra etc. 

Here we focus our attention on Complex Quantum Network Manifolds (CQNMs) of di- 


2 


mension d constructed by gluing together simplices of dimension d. The CQNMs grow 
according to a non-equilibrium dynamics determined by the energies associated to its nodes, 
and have an emergent geometry, i.e. the geometry of the CQNM is not imposed a priori on 
the network manifold, but it is determined by its stochastic dynamics. Following a similar 
procedure as used in several other manuscripts [HHiniiii], one can show that the CQNMs 
characterize the time evolution of the quantum network states. In particular, each network 
evolution can be considered as a possible path over which the path integral characterizing 
the quantum network states can be calculated. Here we show that in d = 2 CQNMs are 
homogeneous and have an exponential degree distribution while the CQNMs are always 
scale-free for d > 2. Therefore for d = 2 the degree distribution of the CQNM has bounded 
fluctuations and is homogeneous while for d = 2 the CQNM has unbounded fluctuations in 
the degree distribution and its structure is dominated by hub nodes. Moreover, in CQNM 
quantum statistics emerges spontaneously from the network dynamics. In fact, here we 
dehne the generalized degrees of the d-faces forming the manifold and we show that the 
average of the generalized degrees of the d-faces with energy e follows different statistics 
(Fermi-Dirac, Boltzmann or Bose-Einstein statistics) depending on the dimensionality 6 of 
the faces and on the dimensionality d of the CQNM. For example in d = 2 the average of 
the generalized degree of the links follows a Fermi-Dirac distribution and the average of the 
generalized degrees of the nodes follows a Boltzmann distribution. In d = 3 the faces of the 
tetrahedra, the links and the nodes have an average of their generalized degree that follows 
respectively the Fermi-Dirac distribution, the Boltzmann distribution and the Bose-Einstein 
distribution. 

Consider a d-dimensional simplicial complex formed by gluing together simplices of di¬ 
mension d, i.e. a triangle for d = 2, a tetrahedron for d = 3 etc. A necessary requirement 
for obtaining a discretization of a manifold is that each simplex of dimension d can be glued 
to another simplex only in such a way that the (d — l)-faces formed by (d — 1)—dimensional 
simplices (links in d = 2, triangles in d = 3, etc.) belong at most to two simplices of 
dimension d. 

Here we indicate with Sd ,5 the set of all d-faces belonging to the d-dimensional manifold 
with d < d. If a (d — l)-face a belongs to two simplices of dimension d we will say that 
it is ’’saturated” and we indicate this by an associated variable with value = 0; if it 
belongs to only one simplicial complex of dimension d we will say that it is ’’unsaturated” 
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and we will indicate this by setting = 1. 

The CQNM is evolving according to a non-eqnilibrinm dynamics described in the follow¬ 
ing. 

To each node i = 1,2... ,iV an energy of the node is assigned from a distribntion 
g{e). The energy of the node is qnenched and does not change dnring the evolntion of the 
network. To every 5-face a G Sd,s we associate an energy e^ given by the snm of the energy 
of the nodes that belong to the face a, 

= ( 1 ) 

zCo: 

At time t = 1 the CQNM is formed by a single d-dimensional simplex. At each time t > 1 
we add a simplex of dimension d to an nnsatnrated {d — l)-face a G Sd^d-i of dimension 
d — 1. We choose this simplex with probability IIq, given by 

Ha = (2) 

where /3 is a parameter of the model called inverse temperature and Z is a normalization 
snm given by 

z= E (3) 

o‘&Sd,d-i 

Having chosen the {d — l)-face a, we glne to it a new d-dimensional simplex containing all 
the nodes of the {d — 1)—face a pins the new node i. It follows that the new node i is linked 
to each node j belonging to a. 

In Fignre [T] we show few steps of the evolntion of a CQNM for the case d = 2, while in 
Fignre 1^ we show examples of CQNM in d = 2 and in d = 3 for different valnes of (3. 

From the dehnition of the non-eqnilibrinm dynamics described above, it is immediate to 
show that the network strnctnre constrncted by this non-eqnilibrinm dynamics is connected 
and is a discrete manifold. 

Since at time t the nnmber of nodes of the network manifold is iV = t + d, the evolntion 
of the network manifold is fnlly determined by the seqnence {ei}i<t+d) and the seqnence 
where e*, for z < d -|- 1 indicates the energy of an initial node, while for i = t' + d 
with F > 1 it indicates the energy of the node added at time t', and where at' indicates the 
(d — l)-face to which the new d-dimensional simplex is added at time t' > 1. 
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Complex Quantum Network Manifold 


The evolution for d=2 



FIG. 1: Few steps of a possible evolution of the CQNM for d = 2. The nodes have different energies 
represented as different colours of the nodes. A link can be saturated (if two triangles are adjacent 
to it) or unsaturated (if only one triangle is incident to each). Starting from a single triangle at 
time t = 1, the CQNM evolves through the addition of new triangles to unsaturated links. 

The dynamics described above is inspired by biological evolntionary dynamics and is 
related to self-organized critical models. In fact the case /? —)■ oo is dictated by an extremal 
dynamics that can be related to invasion percolation [^Hl ITS] , while the case /3 = 0 can be 
identihed as an Eden model BSl on a d-dimensional simplicial complex. 

Here we call these network manifolds Complex Qnantnm Network Manifolds because us¬ 
ing similar arguments already developed in [81- lTO] ITT] it can be shown that they describe 
the evolution of Quantum Network States (see Supplementary Information for details). The 
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quantum network state is an element of an Hilbert space 1-Ltot associated to a simplicial com¬ 
plex of N nodes formed by gluing d-dimensional simplices (see Supplementary Information 
for details). The quantum network state |t/’(^)) evolves through a Markovian non-equilibrium 
dynamics determined by the energies {cj} of the nodes. The quantity Z{t) enforcing the 
normalization of the quantum network state = 1 can be interpreted as a path 

integral over CQNM evolutions determined by the sequences {ei}i<t+d) and {at'}t'<t- In fact 
we have 


= E E W{{ei}i<t+di{oit'}t'<t)i (4) 


where the explicit expression of W{{ei}i<t+d) is given in the Supplementary Infor¬ 

mation. Moreover, Z{t) can be interpreted as the partition function of the statistical me¬ 
chanics problem over the CQNM temporal evolutions. If we identify the sequences {ei}i<t+d, 
and determining Z(t) with the sequences indicating the temporal evolution of the 

CQNM we have that the probability P{{ei}i<t+d, of a given CQNM evolution is 

given by 




^(i) 


(5) 


Therefore each classical evolution of the CQNM up to time t corresponds to one of the paths 
dehning the evolution of the quantum network state up to time t. 

A set of important structural properties of the CQNM are the generalized degrees kd^s 
of its (5-faces. Given a CQNM of dimension d, the generalized degree kd^s^a) of a given 
(5-face (U, (i.e. a G Sd^s) is dehned as the number of (i-dimensional simplices incident to 
it. For example, in a CQNM of dimension d = 2, the generalized degree k 2 ,i{c() is the 
number of triangles incident to a link a while the generalized degree k 2 fi{c() indicates the 
number of triangles incident to a node a. Similarly in a CQNM of dimension d = 3, the 
generalized degrees ^ 3 ^ 2 , ^ 3,1 and k^^ indicate the number of tetrahedra incident respectively 
to a triangular face, a link or a node. If from a CQNM of dimension d one extracts the 
underlying network, the degree Kd{i) of node i is given by the generalized degree kd,o{i) of 
the same node i plus d — 1, i.e. 


Kd{k) = kdfl{i) + d-1. (6) 

We indicate with Pd, 5 {k) the distribution of generalized degrees kd ,5 = k. It follows that 
the degree distribution of the network Pd{K) constructed from the (i-dimensional CQNM is 
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Complex Quantum Network Manifold 
In Dimension d 


cl=2 


d=3 
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FIG. 2: Visualization of CQNM with d = 2 (panel A) and d = 3 (panel B). The colours of the 
nodes indicate their energy {ej} while their size indicates their degree Kd{i)- In d = 2 the degree 
distribution of the CQNMs is a convolution of exponentials, in d = 3 the CQNMs are scale-free and 
the presence of hubs is clearly observable from this visualization. The data shown are for CQNM 
with N = 10^ nodes, /? = 0.2 and Poisson distribution g{e) with average z = 5. 


given by 


Pd(K) = F^,o(^-d+l). (7) 

Let us consider the generalized degree distribution of CQNM in the case /3 = 0. In this case 
the new d-dimensional simplex can be added with equal probability to each unsaturated 
(d — l)-face of the CQNM. Here we show that as long as the dimension d is greater than 
two, i.e. d > 2, the CQNM is a scale-free network. In fact each d-face, with S < d — 1, which 
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has generalized degree kd^si^a) = k, is incident to 


2 + {d - 6 - 2)k 


( 8 ) 


unsaturated {d — l)-faces. Therefore the probability 'n'd, 5 {ci) to attach a new d-dimensional 
simplex to a 5-face a with generalized degree kd^s^a) and with 5 < d — 1, is given by 




2 + {d - 5 - 2)kd,6{,a) 

-f- (d — 5 — 2)kd^&{a'y\ 


(9) 


Therefore, as long as 5 < d — 2, the generalized degree increases dynamically due to an 
effective ’’linear preferential attachment” [ 21 ], according to which the generalized degree of 
a 5-face increases at each time by one, with a probability increasing linearly with the current 
value of its generalized degree. Since the preferential attachment is a well-known mechanism 
for generating scale-free distributions, it follows, by putting 5 = 0, that we expect that as 
long as d > 3 the CQNMs are scale-free. Instead, in the case d = 2, by putting 5 = 0 
it is immediate to see that the probability 7i2fi{c() is independent of the generalized degree 
^ 2 ,o(ct) of the 0 —face (node) a, and therefore there is no ’’effective preferential attachment”. 
We expect therefore nzi that the CQNM in d = 2 has an exponential degree distribution, 
i.e. in d = 2 we expect to observe homogeneous CQNM with bounded fluctuations in the 
degree distribution. These arguments can be made rigorous by solving the master equation 
[T9|, and deriving the exact asymptotic generalized degree distributions for every 5 < d (see 
Supplementary Information for details). For 5 = d — 1 we hnd a bimodal distribution 


^ for k = 1, 


Pd,d-i{k) = 


for k = 2. 

For 5 = d — 2 instead, we hnd an exponential distribution, i.e. 

1 ^ 


Pd,d-2{k) = (^) ^ ^ 1 - 


( 10 ) 


( 11 ) 


Therefore in d = 2, the CQNMs have an exponential degree distribution that can be derived 
from Eq. ( [ll| ) and Eq. 0. Finally for 0<5<d — 2we have the distribution 


PdAk) 


d-l r[l+(d+l)/(d-5-2)] r[k+2/id-5-2)] 

d-5-2 r[l+2/(d-<5-2)] r[A:+l-|-(d+l)/((i-<5-2)] 


, for fc > 1 . 


( 12 ) 


It follows that for 0 < 5 < d — 2 and k S> 1 -|- (d -|- l)/(d — 5 — 2) the generalized degree 
distribution follows a power-law with exponent 7 ^, 5 , i.e. 


PdAk) — Ck for d — 5 > 2, 


( 13 ) 




and 


ld,s = 1 + 


d-l 

d-6-2' 


(14) 


The distribution Pd, 5 {k) given by Eq. (12) is scale-free if an only if 'yd ,5 G (2,3]. Using 


Eq. (14) we observe that for d> 3 and 5 = 0 we observe that the distribution of generalized 


degrees Pd,s{k) is always scale-free. Therefore the degree distribution Pd{K) given by Eq. ([^, 
for large values of the degree K and for d > 3 is scale-free and goes like 


Pd{K) ~ 


(16) 


with 


Id 


2d-3 
d-2' 


(16) 


Therefore, for d = 3 the CQNMs have 73 = 3 and for d —)■ cx) they have power-law exponent 
Id -t 2. 

These theoretical expectations perfectly £t the simulation results of the model as can been 
seen in Figure where the distribution of generalized degrees P 3 ,i{k) and Psflik) observed 
in the simulations for /3 = 0 are compared with the theoretical expectations. 

In the case (3 > 0 the distributions of the generalized degrees depend on the density 
Pd,&{^) of d-dimensional simplices with energy e in a CQNM and are parametrized by self- 
consistent parameters called the chemical potentials, indicated as ^d,& and dehned in the 
Supplementary Information. 

Here we suppose that these chemical potentials iJid,& exist and that the density Pd,s{^) 
is given, and we hnd the self-consistent equations that they need to satisfy at the end of 
the derivation. Using the master equation approach [19] we obtain that for 6 = d — 1 the 
generalized degree follows the distribution 


Pd,d-l{k) — 


while for d = d — 2 it follows 


Pd,d-l (e) 


1 - 


/('-Md.d-i) 


+1 


Pd,d-l(e) ^—Y 


, for /c = 1 
for k = 2, 


(17) 


+1 


Pd,d-2{k) = Y.ePd,d-2{(^) / n k , for P > 1. 


y(e-Md,d-2) 


+ 1 


(18) 
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FIG. 3: The distribution of the (non-trivial) generalized degrees and in dimension 

d = 3 are shown. The star symbols indicate the simulation results while the solid red line indicates 
the theoretical expectations given respectively by Eqs. © and ( |12[ ) . In particular we observe that 
P 3 ,iik) is exponential while P^^^k) is scale-free implying that the CQNM in d = 3 is scale-free. 
The simulation results are shown for a single realization of the CQNM with a total number of 
nodes N = 2 x 10^. 


Finally for 0 < 6 < d — 2 the generalized degree is given by 


PdAk) 


Pd,d-2{e)A^^ 


r[l-|-2/(d-<5-2)-|-exp[/3(e-/jrf_^)]]r[fc-|-2/(d-(5-2)] k > 1 

r[l-|-2/(d-<5-2)]r|^fc-|-l-|-2/(d-5-2)-|-exp[/3(e-/id,j)]] ’ 


(19) 


It follows that also for /3 > 0 the CQNMs in d > 2 are scale-free. Interestingly, we observe 
that the average of the generalized degrees of simplices with energy e follows the Fermi-Dirac 
distribution for 6 = d — 1, the Boltzmann distribution for S = d — 2 and the Bose-Einstein 
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distribution for 6 < d — 2. In fact we have, 


{[kd,d-i - l]\e) = nF{e, Hd,d-i), 

{[kd,d-2 - l]\e) = nz{€,fid,d-2), ( 20 ) 

([fcd,5-l]|e) = AnB{e,Hd,5), for 6 < d - 2. 

where A = {d — 6)/{d — 6 — 2), nz{e, idd,d- 2 ) is proportional to the Boltzmann distribution 
and npie, fid,d-i), nB{e,fid,s) indicate respectively the Fermi-Dirac and Bose-Einstein occu¬ 
pation numbers HU. In particular we have 


nz{e, ^id,d-2) = e Md.d- 2 )^ 

nF{e,fid,d-i) = 

^Bie,f^d,5) = ^ • 


( 21 ) 


These results suggest that the dimension d = 3 of CQNM is the minimal one necessary 
for observing at the same time scale-free CQNMs and the simultaneous emergence of the 
Fermi-Dirac, Boltzmann and Bose-Einstein distributions. In particular in d = 3 the average 
generalized degree of triangles of energy e follows the Fermi-Dirac distribution, the average 
of the generalized degree of links of energy e follows the Boltzmann distribution, while 
the generalized degree of nodes of energy e follows the Bose-Einstein distribution. Finally 
the chemical potentials fid, 5 , if they exist, can be found self-consistently by imposing the 
condition 


lim TsiIm = 

t^oo Nd,S d -I- 1 


( 22 ) 


dictated by the geometry of the CQNM, which implies the following self-consistent relations 
for the chemical potentials fid ,5 


J depd,d-i{e)nF{e, fdd,d-i) — 

J depd,d-i{e)nz{e, fid,d- 2 ) = 
j depd,s{e)nBie, pd,s) = ^ ^ 


for 


6 <d- 2 . 


(23) 

(24) 

(25) 


In Figure we compare the simulation results with the theoretical predictions given by 


Eqs. (20) Ending very good agreement for sufficiently low values of the inverse temperature f3. 
The disagreement occurring at large value of the inverse temperature fd is due to the fact that 
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FIG. 4: The average of the generalized degrees of d-faces of energy e, in a CQNM of dimension d = 3, 
follow the Fermi-Dirac nF{e, ^ 3 ^ 2 ), the Boltzmann ras(e,/i 3 ^i) or the Bose-Einstein distribution 
according to Eqs. — (25) as long as the chemical potential ^d,s is well-defined, i.e. 


for sufficiently low value of the inverse temperature /3. Here we compare simulation results over 
CQNM oi N = 10^ nodes in d = 3 and theoretical results for /3 = 0.01,0.1,1. The CQNMs in 
the figure have a Poisson energy distribution g^e) with average z = 5. The simulation results are 
averaged over 50 CQNM realizations. 


the self-consistent Eqs. (23) — (25) do not always give a solution for the chemical potentials 


fid, 5 - In particular the CQNM with d > 3 can undergo a Bose-Einstein condensation when 


Eq. (25) cannot be satished. When the transition occurs for the generalized degree with 


(5 = 0, the maximal degree in the network increases linearly in time similarly to the scenario 
described in [27] . 

In summary, we have shown that Complex Quantum Manifolds in dimension d > 2 are 
scale-free, i. e. they are characterized by large fluctuations of the degrees of the nodes. 
Moreover the d-faces with 6 < d follow the Fermi-Dirac, Boltzmann or Bose-Einstein distri¬ 
butions depending on the dimensions d and 6. In particular for d = 3, we hnd that triangular 
faces follow the Fermi-Dirac distribution, links follow the Boltzmann distribution and nodes 
follow the Bose-Einstein distribution. Interestingly, we observe that the dimension d = 3 
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is not only the minimal dimension for having a scale-free CQNM, but it is also the mini¬ 
mal dimension for observing the simultaneous emergence of the Fermi-Dirac, Boltzmann or 
Bose-Einstein distributions in CQNMs. 
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SUPPLEMENTARY INFORMATION 


INTRODUCTION 

In this supplementary information we give the details of the derivation discussed in the 
main text. In Sec. II we dehne Complex Quantum Network Manifolds (CQNMs); in Sec. Ill 
we discuss the relation between the CQNM and the evolution of quantum network states; 
hnally in Sec. IV we dehne the generalized degrees, and we derive the generalized degree 
distribution in the case /d = 0 and /d > 0. 

COMPLEX QUANTUM NETWORK MANIFOLDS 

Here we present the non-equilibrium dynamics of Complex Quantum Network Manifolds 
(CQNMs). This dynamics is inspired by biological evolution and self-organized models and 
generates discrete manifolds formed by simplicial complexes of dimension d. In particular 
CQNMs are formed by gluing d-simplices along (d— l)-faces, in order that each (d— l)-face 
belongs at most to two d-dimensional simplices. 

Let us indicate with Sd,s the set of all d-faces with 6 < d belonging to the d-dimensional 
CQNMs. A (d— l)-face a G Sd,d-i is ’’saturated” if it belongs to two simplices of dimension 
d, whereas it is ’’unsaturated” if it belongs only to a single d-dimensional simplex. We 
will assign a variable = 0,1 to each face a G Sd^d-i, indicating either that the face is 
unsaturated (.^q, = 1) or that the face is saturated (^q, = 0). 

Moreover, to each node i we assign an energy drawn from a distribution g{e) and 
quenched during the evolution of the network. To every d-face a G Sd,s we associate an 
energy given by the sum of the energy of the nodes that belong to a, 

ea = '^ei. (S-1) 

iGa 

At time t = 1 the CQNM is formed by a single d-dimensional simplex. At each time t > 1 
we add a simplex of dimension d to an unsaturated (d — l)-face a G Sd,d-i chosen with 
probability Hq, given by 

n„ = (S-2) 
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where /3 is a parameter of the model called inverse temperature and Z is a normalization 
snm given by 

Z= (S-3) 

a&Sd,d-l 

Having chosen the {d — l)-face a, we glne to it a new d-dimensional complex containing all 
the nodes of the face a pins the new node i. It follows that the new node i is linked to each 
node j belonging to a. 

Since at time t = 1 the nnmber of nodes in the CQNM is A^(l) = d + 1, and at each time we 
add a new additional node, the total nnmber of nodes is A^(t) = t-j-d. The CQNM evolntion 
np to time t is fnlly determined by the seqnences {ei}i<t+d, where et'+d indicates 

the energy of the node added to the CQNM at time t' > 1, with i < d + 1 indicates the 
energy of an initial node i of the CQNM, and at/ indicates the {d — l)-face to which the new 
d-dimensional complex is added at time t'. A similar dynamics for simplicial complexes of 
dimension d = 2 has been proposed in |40j . 


QUANTUM NETWORK STATES 


The network Hilbert space 


Following an approach similar to the one nsed in ” Qnantnm Graphity” and related models 
[HHiniil], in this section we associate an Hilbert space 'Htot to a simplicial complex of N 
nodes formed by glning together d-dimensional simplices along (d — l)-faces. The Hilbert 
space Titot is given by 

N p p 

'Htot = Hnode Hd,d-1 'Hd,d-1, (S-4) 


with P 



indicating the maximnm nnmber of (d — l)-faces in a network of N nodes. 


Here an Hilbert space Hnode is associated to each possible node i of the simplicial complex, 
and two Hilbert spaces Hd,d-i and Hd,d-i are associated to each possible (d — 1)—face of 
a network of N nodes. The Hilbert space Hnode is the one of a fermionic oscillator of 
energy e, with basis {|oi,e)}, with ot = 0,1. We indicate with h\{e),hi{e) respectively the 
fermionic creation and annihilation operators acting on this space. The Hilbert space Hd,d-i 
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associated to a (d— l)-face a is the Hilbert space of a fermionic oscillator with basis {|aQ,)}, 
with Oq, = 0,1. We indicate with cj^jCa respectively the fermionic creation and annihilation 
operators acting on this space. Finally the Hilbert space 'Hd,d-i associated to a (d — l)-face a 
is the Hilbert space of a fermionic oscillator with basis {|riQ,)}, with Ua = 0,1. We indicate 
with hjj, ha respectively the fermionic creation and annihilation operators acting on this 
space. 

A quantum network state can therefore be decomposed as 

IV’(^)) = C'({0i,ei,a„,n„}) JJ |0i,ei) JJ M \na), (S-5) 

{oi,ei,aa,na} i a^Qd,d-i(N) 

where with Qd,d-i{N) we indicate all the possible (d — l)-faces of a network of N nodes. 

The node states |oi, e) are mapped respectively to the presence (|oj = 1, e)) or the absence 
(|oj = 0,e)) of a node i of energy e* = e in the simplicial complex. The quantum state 
\aa = 1) is mapped to the presence of the (d — 1)—face a G Sd,d-i in the network while the 
quantum state Ioq = 0) is mapped to the absence of such a face. Moreover, when Oq, = 1, 
the quantum number = 1 is mapped to a saturated (d — 1)—face a, i.e. a is incident 
to two d-dimensional simplices, while the quantum number Uq = 0 is mapped either to an 
unsaturated (d — 1)—face a (if also Oq, = 1) or to the absence of such a face (if Oq, = 0). 

Markovian evolution of the quantum network states 

As already proposed in the literature mm, here we assume that the quantum network 
state follows a Markovian evolution. In particular we assume that at time t = 1 the state is 
given by 

.. d+l 

iV'(i)) = —n 

where 2{1) is hxed by the normalization condition ('0(1)|'0(1)) = 1. The quantum network 
state is updated at each time f > 1 according to the unitary transformation 

m)) = Utmt-i)) (s-7) 
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with the unitary operator Ut given by 


Ut — \ l ^ ^ ^ ' 

Q:£Qd,d—1 




n 

a' (t+d,a) 


J 


^i4ca(S-8) 


where J^(i, a) indicates the set of all the {d — l)-faces a' formed by the node i and a subset 
of the nodes in a G Qd,d-i{N) and Z{t) is hxed by the normalization condition 


= 1- 


(S-9) 


Path integral characterizing the quantum network state evolution 


The quantity Z(t) is a path integral over CQNM evolutions determined by the sequences 


{ei}i<t+d, In fact, using the normalization condition in Eq. (S' — 9) and the evolu¬ 


tion of the quantum network state given by Eqs. (S — 7), (S — 8) we get 


W{{ei}i<t+d, 


(S-10) 


where W{{ei}i<t+d, is given by 

t+d 


Z=1 


t'<t 


where the terms aa{t) and riaii) that appear in Eq. (S — 11) can be expressed in terms of 
the history as 

= X] X] S[a,d] + ^ 6[a,d], 


n, 


t'=2 \a&T{a^,,t'+d) 
t-1 

,{t) = '^5[at',a]. 


asQd_d_i((i+i) 


(S-12) 


t'=2 


We note that Z{t) can also be interpreted as the partition function of the statistical 
mechanics problem over possible evolutions of CQNM. In fact, the CQNM evolution is 
determined by the sequences {ati}t'<t where ef+d indicates the energy of the node 

added at time f' > 1 and e* for i < d -|- 1 indicates the energy of a node i in the CQNM 
at time f = 1, and at> indicates the (d — l)-face to which we attach the new d-dimensional 
simplex at time t'. The probability P{{ei}i<t+d) of a given evolution is given by 

W({ 




m 


(S-13) 
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where W{{ei}i<t+d, is given by Eq. (S' — 11) and Z{t) is fixed by the condition Eq. 


(S' — 9). This implies that the set of all classical evolutions of the CQNM fully determine the 
properties of the quantum network state evolving through the Markovian dynamics given 


byEq. (^-7). 


GENERALIZED DEGREES 

Definition of generalized degrees 

A set of important structural properties of the CQNM are the generalized degrees fed ,5 (a) 
of the (5-faces in a (i-dimensional CQNM. Given a CQNM of dimension d, the generalized 
degree kd^s{a) of a given ^-face a, (i.e. a G Sd^s) is defined as the number of (i-dimensional 
simplices incident to it. If we consider the adjacency tensor a of elements Oq = 1 if the 
d-dimensional complex is part of the CQNM and otherwise zero, Oq, = 0, the generalized 
degree of a h-face a' is given by 

kd,5{a') = ^ tta- (S-14) 

cx\a'Ga 

For example, in a CQNM of dimension d = 2, the generalized degree ^ 2 , 1 ( 0 ;) is the number 
of triangles incident to a link a while the generalized degree ^ 2 , 0 ( 0 ^) indicates the number 
of triangles incident to a node a. Similarly in a CQNM of dimension d = 3, the general¬ 
ized degrees ^ 3 ^ 2 , ^ 3,1 and k^ Q indicate the number of tetrahedra incident respectively to a 
triangular face, a link or a node. 


Distribution of Generalized Degrees for /3 = 0 

Let us define the probability 7rd,5{ci) that a new (i-dimensional simplex is attached to a ^- 
face a. Since each (i-dimensional simplex is attached to a random unsaturated {d— 1)—face, 
and the number of such faces is {d — l)t, we have that for 6 = d — 1 

f Tj^TT for k = 1, 

0 for fc = 2. 
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Let us now observe that each 5-face, with 5 < d—1, which has generalized degree kd, 5 {a) = k, 
is incident to 


2 + {d-5-2)k, 


(S-15) 


unsaturated (d — l)-faces. 

In fact, is is easy to check that a 5—face with generalized degree k^^s = /c = 1 is incident to 
d—d unsaturated (d—l)-faces. Moreover, at each time we add to a 5-face a new d-dimensional 
simplex, a number d — 5 — 1 of unsaturated (d — 1)—faces are added to the 5-face while a 
previously unsaturated (d — l)-face incident to it becomes saturated. Therefore the number 
of (d — l)-unsaturated faces incident to a 5-face of generalized degree kd,s = k follows Eq. 


{S — 15). We have therefore that the probability Tid^s^ci) to attach a new d-dimensional 
simplex to a 5-face a with 5 < d — 1 and generalized degree kd^s{(^) is given by 

2 + {d - 5 - 2)kd,5{a) 




P + {d - S - 2)ka.s(a')y 


(S-16) 


where for large times t S> 1, 


^ [2 + (d-5-2)A;rf,5(a')] =2 I I f+(d-5-2) I I f= (d-1) ( \ t. 


a'eSd,5 


5 + 1 


5 + 1 


From Eq. (5' — 16) if follows that, as long as 5 < d — 2, the generalized degree follows a 


’’preferential attachment” mechanism 

Moreover, the average number Ud^sik) of 5-faces of generalized degree kd^s = k that 
increases their generalized degree by one at a generic time t > 1, is given by 




(S-17) 


d 


where mds+i = I indicates the total number of 5—faces incident to the d- 

5 + 1 

dimensional simplex added at time t. This implies that for 5 = d — 1, and large times 
f » 1, nd, 5 ik) is given by 


nd,d-i{k) = 


(d — l)t 


dfc,i 


where 6x^y indicates the Kronecker delta, while for 5 < d — 1 it is given by 


nd,s(,k) = 


2 + {d-6-2)k 
(d — l)t 


(S-18) 


(S-19) 


21 







Using Eqs. (S' — 18) — (S' — 19) and the master equation approach [19] , it is possible to derive 
the exact distribution for the generalized degrees. We indicate with N\^{k) the average 
number of h-faces that at time t have generalized degree = k during the temporal 
evolution of a d-dimensional CQNM. The master equation [19] for N^g{k) reads 


- l)Nls{k - 1)(1 - 4,i) - n^AANlsik) + md,sSk,i (S-20) 

d 

with k > 1. Here 5x,y indicates the Kronecker delta, and rud^s = j I is the number of 

h-faces added at each time t to the CQNM. The master equation is solved by observing that 
for large times t1 we have N\^^{k) ~ md^stPd, 5 {k) where Pd^sik) is the generalized degree 
distribution. For d = d — 1 we obtain the bimodal distribution 


for A; = 1 


Pd,d-l{k) — 


j for fc = 2 

For 6 = d — 2 instead, we hnd an exponential distribution, i.e. 

Pd,d- 2 {k) = for fc > 1 . 

Finally for 0 < 6 < d — 2 we have the distribution 

p ( _ d-l r[l+(d+l)/(d—2)] V[k+2/{d—5—2)] r > 1 

d-5-2 r[l+2/(d-<5-2)] r[A;+l+(d+l)/((i-<5-2)] ’ ft. ^ l. 


(S-21) 


(S-22) 


From Fq. (S — 22) it follows that for 0 < d < d — 2 and /c S> 1 the generalized degree 


distribution follows a power-law with exponent 'jd^s, i.e. 

Pd,s{k) ~ Ck~'^^’^ for d — d > 2, 


(S-23) 


and 


Id,5 = 1 + 


d — 1 
d-5-2' 


(S-24) 


Therefore the generalized degree distribution Pd, 5 {k) given by Fq. (S' —22) is scale-free, 
i.e. it has diverging second moment {{kd, 5 )‘^), as long as G (2,3]. This implies that the 
generalized degree distribution is scale-free for 

d-3 


d < 


(S-25) 
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Distribution of Generalized Degrees for /3 > 0 


In the case /? > 0 the distribution of the generalized degrees Pd, 5 {k) are convolutions 
of the conditional probabilities Pd, 5 {k\e) that 5-faces with energy e have given generalized 
degree kd,s = k. Here we derive the distribution of generalized degrees Pd, 5 {k) for different 
values of 5 and d as a function of the inverse temperature /3. The procedure for hnding these 
distributions is similar for every value of 5. 

First we will determine the master equations [19] for the average number iVj^(fc|e) of 5- 
faces of energy e that have generalized degree k at time t. Then we will solve these equa¬ 
tions, imposing the scaling, valid in general for growing networks, given by Nd^ 5 {k\e) = 
'^d, 5 Pd, 5 {^)tPd, 5 {k\e), where Pd,s{^) is the probability that a 5-face has energy e and md,s are 
the number of 5-faces added at each time to the CQNM. The master equations will also 
depend on self-consistent parameters pd,s called chemical potentials that need to satisfy self- 
consistent equations for the derivation to hold. 

Let us consider hrst the case 5 = 5—1. The average number N\^^{k\e) of 5 = (d— l)-faces of 
energy e that at time t have generalized degrees kd^s = k follows the master equation given 
by 

= 2|t) - NX^_dk = 2|t) = = l|t), 


^ l\0 - Ni,_dk ^ l\0 


--^K.d-i(k = 1|€) + (S-26) 


where pd 4 -i{e) is the probability that a (d — l)-face added to the network at a generic 


time t 3> 1 has energy e and rud^d-i 



d is the number of 5-faces added to 


the network at each time t. In order to solve this master equation we assume that the 
normalization constant Z <x.t and we put 


g-/3/2d,d-i = lim-. (S-27) 

t^oo t 

This is a self-consistent assumption that must be verihed by the solution of Eqs. ( [S' — 26 ). 
Moreover we observe that at large times Nj^d-iik\^) — ’md,d-iPd,d-i{^)Pd,d-i{k\e). Here 
Pd^d-i{k\e) indicates the asymptotic probability that a (d — l)-face a with energy e has 
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kd,d-i{(^) = k. With these assumptions, we can solve Eqs. (S' — 26) Ending 


Pd4-iik = l\e) = Pd,d-i(e) 


e/3h-Md,d-i) + 1 - Pd4-i{^)[^ - nF{e, ^id4-i)] 


Pd 4 -i{k — 2\e) — — Pd,d-i(e)ni;'(e,/irf,rf_i), 


(S-28) 


where np{e, fid,d-i) is the Fermi-Dirac occupation number with chemical potential fid,d-i, 

i.e. 

1 


npi^, Pd,d-i) — 


g/3(e 


(S-29) 


Using Eqs. (S' — 28), and performing the average (fcd,d-i|e) of the generalized degree kd^-i 
over {d — l)-faces of energy e, one can easily hnd that 

1 


{[W-i - Ilk) = ^ 1 


(S-30) 


This result shows that the average of generalized degrees of (d — l)-faces of energy e is 
determined by the Fermi-Dirac statistics with chemical potential ^d,d-i- 
Let us now consider the case 6 = d — 2. In this case we assume that asymptotically in time 
we can dehne the chemical potential fid,d -2 as 


g/^Md,d — 2 — —1 


t^OO 


SaeQd.d_2(t) Sa'eQd.d_i(t)|aCa' ^ °‘\a'd{kd4-2{oi) ^ k) 

XlaeQd,d_2W k) 


(S-31) 


In this assumption, the master equations [H] for the average number N\^_ 2 {k\e) of (d — 2)- 
faces with energy e and generalized degree fc > 1, read 


Nl^d-2{k\e) - Nl,_2{k\e) = 


p — 2 ) p — 2 ) 

■" "■ ‘ ^— N‘(m 


Wi,_2(;;-lk)|l-4,i] 


+''n^d,d-2Pd,d-2{.^)5k,l^ 


(S-32) 


where md 4-2 = d(d — l)/2 is the number of (d — 2)—faces added at each time t to 
the CQNM, Pd,d- 2 {.^) is the probability that such faces have energy e, and indicates 
the Kronecker delta. In the large network limit t S> 1 we observe that Nj^d- 2 ik\^) — 
tmd,d- 2 Pd,d- 2 {^)Pd,d- 2 {k\e) where Pd,d- 2 {k\e) indicates the probability that a (d —2)—face of 


energy e has generalized degree k. Solving Eq. (S' — 32) we get. 


Pd,d- 2 {k\e) = Pd,d- 2 (e) 


ojdp — 2 ) 


|^e/3(e-Md,d-2) _|_ xj 


k ’ 


(S-33) 
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for A; > 1. Therefore, summing over all the values of the energy of the nodes e we get the 
full degree distribution P{k) 


Pd,d-2{k) — Pd,d-2(e) 


3/3(e—A‘d,d-2) 




k ’ 


(S-34) 


for k > 1. Using Eqs. (S' —34), and performing the average {kd,d- 2 \^) of the generalized 
degree kd 4-2 over (d — 2)-faces of energy e, one can easily hnd that 


{[kd,d-2 - l]|e) = e 


^z(e, fJ'd,d-2), 


(S-35) 


where we have indicated with nzi^, ^d,d- 2 ) the Boltzmann distribution 


l^d,d-2) 


= g-/3(e-Md.d-2) 


(S-36) 


This result shows that the average of generalized degrees of (d — 2)-faces of energy e is 
determined by the Boltzmann statistics with chemical potential iid,d- 2 - 
Let us hnally consider the case 6 < d — 2. In this case we assume that asymptotically in 
time we can dehne the chemical potential fid,s given by 




t^OO 


Assuming that the chemical potential ^d,s exists, the master equations [19] for the average 
number A"^^(A;|e) of d-faces with energy e and generalized degree A; > 1 read 

ivijtfck) - ivtj(ik) = -1 + Vid - s - 


e-/3h-Md.b[fc + 2/(d-d-2)] 


Nisik\e) 


+'md,sPd,s{^)dk,i, 


(S-38) 


where md,s = 



is the number of d—faces added at each time t to the CQNM, Pd,s{,^) is 


the probability that such faces have energy e, and 6x,y indicates the Kronecker delta. In the 
large network limit A S> 1 we observe that Nlj^§{k\e) ~ tmd^ 5 Pd, 5 {^)Pd, 5 {k\e),where Pk^s{k\e) is 


the probability that a d—face of energy e has generalized degree k. Solving Eq. (S' — 38) we 
get. 


PdAk) 


EeP<i,d-2(e)e^(" 


r[l+2/(d-5-2)+exp[/3(6-Md,d)]]r[fc+2/(d-,5-2)] 
r[l+2/(d-5-2)]r[fc+l+2/(d-5-2)+exp[/3(e-/2d,5)]] ’ “ ' ^ ’ 
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Using Eqs. (S' — 39), and performing the average {kd,i\e) of the generalized degree over 
5-faces of energy e, one can easily hnd that 


{[kd,5-l]\e) = = Annie, Hd, 5 ), (S-40) 

where A = {d — 6)/{d — 6 — 2) and where ^^(e, fid,s) indicates the Bose-Einstein occupation 
number with chemical potential fid,s, he. 

= (s-41) 

This result shows that the average of generalized degrees of 5-faces of energy e is determined 
by the Bose-Einstein statistics with chemical potential fid, 5 - Finally the chemical potentials 
^d, 5 , if they exist, can be found self-consistently by imposing the condition 


lim 

t^OO 


Nd,5 



d -|- 1 

5n’ 


(S-42) 


dictated by the geometry of the CQNM. In fact at each time t we add to the network 
d 


md,5 = 


new 5-faces and we increases the sum of the generalized degree kd^s by the 


amount 


d -|- 1 

5 + 1 

for the chemical potentials fid, 5 , 



. Imposing Eq. (S' — 42) implies the following normalization constraints 


/ depd,d-i{e)nF{e, fJ^d,d-i) = 

/ depd,d-i{€)nz{e, fid,d-2) = 

/ depd,&{e)nB{e, pd,5) = for 5 < d - 2. 


(S-43) 


<5+1 


For small values of /3, these equations have a solution that converges for /3 —>■ 0 to the /? = 0 
solution discussed in the previous subsection. As the value of fd increases it is possible that 
the chemical potentials pd,s become ill-dehned and do not exist. In this case different phase 
transitions can occur. For the case d = 2 these transitions have been discussed in detail in 
HU. For d > 2 we observe that the network might undergo a Bose-Einstein phase transition 


for values of the inverse temperature for which Eq. (5' — 43) cannot be solved in order to 
hnd the chemical potential pd,s- The detailed discussion of the possible phase transitions in 
CQNM is beyond the scope of this work and will be the subject of a separate publication. 
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Mean-field treatment of the case j] > 0 


It is interesting to characterize the evolution in time of the generalized degrees using 
the mean-held approach [12]. This approach reveals other aspects of the model that are 
responsible for the emergence of the statistics determining the distribution of the generalized 
degrees. Let us consider separately the mean-held equations determining the evolution of 
the generalized degrees of 5—faces with 6 = d— l,d — 2 or with 6 < d — 2. 

The (d — l)-faces can have generalized degree kd^d-i that can take only two values = 

1, 2. The indicator of a (d — 1)—face a with generalized degree = k is given by 

^a = 2-kd,d-i{a). (S-44) 


In fact for k = 2 the face is saturated and ■Ca = 0 while for k = 1 the face is unsaturated, 
therefore = 1- The mean-held approach consists in neglecting huctuations, and identifying 
the variable (evaluated at time t for a d—face a arrived in the CQNM at time ta) with 
its average t®) over all the CQNM realizations. The mean-held equation for is 

given by 


dja _ e 
dt Z 


(S-45) 


with initial condition ^a{ta, ta) = 1 where ta is the time at which the (d —1)—face is added to 


the CQNM. The dynamical Eq. (S' — 45) is derived from the dynamical rules of the CQNM 


evolution. In fact, at each time one (d — 1)—face is chosen with probability IIq, given by 


Eq. (S' — 2). This face becomes unsaturated and glued to the new d-dimensional simplex. 


Therefore at each time ^a indicating the average of ^a decreases in time by an amount given 
by IIq,. Assuming that for large time t we have Z ~ where the chemical potential 


hrf.d-i is dehned in Eq. (S' — 27), it follows that the solution of the mean-held Eq. (S' — 45) 
is given by 


^a{,t,ta) = ( J 


The average of over all (d — 1)—faces a with energy e, i.e. (CqIc) is given by 


(S-46) 


^qIc) = {[2 - kd4-i]\e) = - dtai — 


t 


tr 


t 


-r + o(f-Q. (S-47) 

1 + e/3(e-/dd,d-l) ' ^ ^ ^ ^ 
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Therefore we obtain also in the mean field approximation, that the average generalized 
degree of (d — 1)—faces with energy e, for f S> 1 satishes, 




l]k) = 




(S-48) 


_j_ 

where npi^, fid,d-i) indicates the Fermi-Dirac occupation number with chemical potential 
Let us consider now the generalized degree of {d — 2)-faces using the mean-held approxima¬ 


tion. Assuming that the chemical potential ^id,d -2 dehned in Eq. (F — 31) is well dehned, it 


is possible to write down the mean held equation for the average kd,d- 2 {oi) of the generalized 
degree of the (d — 2)-face a over CQNM realizations. The solution of this equation will 
provide the evolution of the average of the generalized degree of a (d — 2)— face a with 
energy at time t, given that the face a is added to the CQNM at time ta, i.e. the solution 
will specify the function kd,d- 2 {,C() = kd,d- 2 it,ta)- The mean-held equation is given by 

dkd4-2{oi) 


dt 


(S-49) 


with solution ^^,^- 2 ( 0 ) = kd,d- 2 it,ta) and initial condition kd,d- 2 ita,ta) = 1. It follows that 
kd,d- 2 (t,ta) evolves in time as 


kd,d- 2 {t,to,) = e X. 


(S-50) 


The average {kd,d -2 — l|e) of the generalized degree minus one of (d — 2)—faces with energy 
e is given by 

{[kd,d-2 - l]|e) = ([kd,d-2 - l]|e^ + Oi\\\t/t). 

Therefore we obtain also in the mean held approximation, that the average generalized 
degree of (d — 2)—faces with energy e, for t S> 1 satishes. 


{[kd,d-2 - l]|e) = e' 


- —2) - 


— /id,d-2), 


(S-51) 


where nz{e^ iJLd 4 - 2 ) is proportional to the Boltzman distribution at temperature T = 1/fd. 
Finally the mean-held equation for the generalized degrees of d-faces with d < d —2, is given 
by 

dkd,&{a) [fcd,5(a) + 2/(d - d - 2)] 


dt 


(S-52) 








where one considers the solution kd^sioi) = kd, 5 {t,ta) with initial condition kd^&{ta,ta) = 1. 
Here kd^s indicates the average over CQNM realizations of the average degree of the 5—face 


a, and ^d ,5 is the chemical potential dehned in Eq. (S' — 37). The mean-held solution of Eq. 


(S — 52) is given by 


kd,S (^) to) 


d-6-2 


= A - 


„-/3(eQ:-Md,5) 


(S-53) 


where A = 1 + 2/{d — 5 — 2) = [d — 5)/{d — 5 — 2). The average {[kd^d -2 — l]|e) of the 
generalized degree minus one of the (d — 2)—faces with energy e is given by 


{[kd,5 - 1]|6) = {[kd,5 - l]|e) = H- 



f t\ 

/ dt(y^ 

^ -1 

Jl 

\taj 


= A 




+ 0{t 




Therefore we obtain also in the mean held approximation, that the average generalized 
degree of 5—faces with energy e, and (5<d — 2, fort3>l satishes, 

1 


{[kd,s - l]|e) = A- 


= Ansie, iJ,d,5), 


(S-54) 


g/3(e t^d,5) — 

where nB{e,fid, 5 ) is proportional to the Bose-Einstein distribution at temperature T = 1//9. 


As a hnal remark we note that the mean-held Eqs. (S — 45), (S — 49), (S — 52) can be 
written by the single equation 




(S-55) 


with tti = —1, 0,1. 


In fact Eq. (S — 45) is equal to Eq. (S — 55) where ai = —1 and y = ^a- Eq. (S — 49) 


is equal to Eq. (S — 55) with a* = 0 and y = kd^-i- Finally Eq. (S — 52) is equal to Eq. 
(S — 49) with Oj = 1 and y = kd^s + 2/{d — 5 — 2). 
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